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ABSTRACT 

Two different mechanisms may act to induce quasi-periodic pulsations (QPP) in whole- 
disk observations of stellar flares. One mechanism may be magneto-hydromagnetic 
(MHD) forces and other processes acting on flare loops as seen in the Sun. The other 
mechanism may be forced local acoustic oscillations due to the high-energy particle 
impulse generated by the flare (known as ‘sunquakes’ in the Sun). We analyze short- 
cadence Kepler data, of 257 flares in 75 stars to search for QPP in the flare decay branch 
or post-flare oscillations which may be attributed to either of these two mechanisms. 
About 18 percent of stellar flares show a distinct bump in the flare decay branch of 
unknown origin. The bump does not seem to be a highly-damped global oscillation 
because the periods of the bumps derived from wavelet analysis do not correlate with 
any stellar parameter. We detected damped oscillations covering several cycles (QPP), 
in seven flares on five stars. The periods of these oscillations also do not correlate with 
any stellar parameter, suggesting that these may be a due to flare loop oscillations. We 
searched for forced global oscillations which might result after a strong flare. To this 
end, we investigated the behaviour of the amplitudes of solar-like oscillations in eight 
stars before and after a flare. However, no clear amplitude change could be detected. 
We also analyzed the amplitudes of the self-excited pulsations in two S Scuti stars and 
one 7 Doradus star before and after a flare. Again, no clear amplitude changes were 
found. Our conclusions are that a new process needs to be found to explain the high 
incidence of bumps in stellar flare light curves, that flare loop oscillations may have 
been detected in a few stars and that no conclusive evidence exists as yet for flare 
induced global acoustic oscillations (starquakes). 
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1 INTRODUCTION 


Photometric observations from the Kepler spacecraft have 
revealed flares in a considerable number of stars. These oc¬ 
cur not only in cool M dwarfs, but in K, G, F and A stars as 
well. The numbers of B stars observed by Kepler are too 
few to allow detection of flares. Because even the weak¬ 
est visible stellar flares are thousands of times more ener¬ 
getic t han a large solar flare, t hese are often called “super¬ 
flares” . IWalkowic.z et all (1201 lT ) were the first to identify 373 
flare stars out of ~ 230 00 coo l dwa rfs in the Kepler field. 
Subsequently, iMaehara et al.l d2012lf disc overed 148 solar- 
type s tars with superflares. More recently, Shibavam a et all 
(|2013|) found superflares on 279 G dwarfs. fealonal ( 20121 ') and 
iBalonal d2013l i discovered several A stars with superflares 


which cannot be attributed to flares in a cool companion. It 
appears that A stars may have spots and flares in spite of 
the lack of significant convection. In fact, the incidence of 
flares on A st ars is not much lower than in G and F stars 
dBalonall2015h . In these hot stars the magnetic field may be 
a res ult of the Tayler instability in a different ially-rotating 
star dSpruitll2002l :' lMullan fe Mac.Donaidll2005i l. 

Whereas solar flares emit predominantly in several dis¬ 
crete chromospheric and UV lines of highly-ionized elements, 
the optical flux in stellar flares appears to be distributed 
much like the continuum of an A or B star dKowalski et aid 
l2013ll . White-light emission in solar flares is difficult to de¬ 
tect, though it is possible that it mig ht contribute a larg e 
fraction of the total flare radiation dKretzschman 12011( 1. 
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Even the strongest solar flares are barely detected in space 
observations of total solar irradiance and the Sun would not 
be detected as a flare star by Kepler. It is thought that solar 
flares arise from the energy released by magnetic reconnec¬ 
tion. Although superflares are typically 10 I * * * * 6 times more ener¬ 
getic than large solar flares, it is possi ble that the magn etic 
reconnection model may still apply dShibata et al.ll2013n . 

Nearly all Kepler photometry has been obtained with 
30-rnin long-cadence (LC) exposures, which means that only 
flares which have a long duration can be detected. A few 
thousand stars were observed with 1-min short-cadence (SC) 
exposures, though only for a relatively short time. Whereas 
LC data are nearly continuous over the four-year time span 
of Kepler observations, the SC data cover typically only a 
few months for any star. The advantage of S C data is that it 
allows flares of short duration to be detected. iBalonal (|2015l f 
discovered 3140 flares in 209 stars of all spectral types ob¬ 
served in SC mode. In this paper the flare in a particular star 
is identified by its sequence number. For example 002300039- 
028 is the 28-th flare in KIC 2300039 listed in the catalogue 
of IBalonal (l2015l ~l. 


We know that quasi-periodic p ulsations (QPP) of coro¬ 
nal loops occur in some solar f lares dNakariakov fe Melnikovl 
120091 : iNakariakov et ahl l20ld ~l. The mechanism giving rise 
to QPP is not fully understood at prese nt. One possibility 
pointed out by [M cLau ghl in et al. l (120121 ) is that of oscilla¬ 
tory reconnection. The behaviour of oscillatory reconnection 
is similar to a damped harmonic oscillator and may play a 
role in generating QPP. 

Analysis of QPP can provide information on the 
flare coronal environmen t and magnetic field stre ngth by 
the use of seismology (iKuprivanova et al.l l2013l 'l. How¬ 
ever, solar flare QPP are only observed in Ha, extreme 
UV lines and in radio and X-ray and gamma-ray ob¬ 
servations. There are no recorded observations of QPP 
in white-light solar flares. T he reported optical observa¬ 
tions of QPP in stellar flar es (|Rodono||l974|; Zhilyaevet_alJ 
2000L Mathloudakis_et_aL|_ 200fi Contadakis et alT ~ 20ld : 


Qian et al.l |2012| : lAnfinogentov et al .1120131 ). which have pe¬ 
riods ranging from a few seconds to tens of minutes, are 
whole-disk essentially white-light observations. They can¬ 
not therefore be directly compared with solar flare QPP. It 
is also far from certain if white-light stellar flares are gener¬ 
ated in the same way as solar flares. 

Solar and stellar QPP have previously been used to pro¬ 
vide estimates for several flare parameters. An example of 
QPP with large amplitude and duration in a solar flare ob- 
served in X - ray an d microwave radi o bursts in described by 

I Kane et al.l (Il983l l. More recently, IVan Doorsselaere et al.l 
( 201 ill used X-ray observations of two oscillation modes in a 

single solar flare to estimate the plasma-beta and the density 

contrast of the flaring loop. The wave m ode number was also 

estima ted from the observed periods. lAnfinogentov et al.l 
(l2013ll analyzed the oscillatory signal in the decay phase of 

the U-band light curve of a flare in the dM4.5e star YZ CMi. 
The observational signature is typical of the longitudinal os¬ 
cillations observed in solar flares at extreme ultraviolet and 

radio wavelengths and is associated with standing slow mag¬ 
netoacoustic waves. They therefore suggest that the QPP 
in this stellar superflare may be of a similar nature to so¬ 
lar QPP. A well pronounced QPP has been reported dur¬ 
ing an very energetic flare on the RS CVn binary II Peg 


dMathioudakis et al.ll2003h . The QPP has a long period lead¬ 
ing to the peak. 

The dynamic impact in the photosphere caused by 
a solar flare is called a “sunquake”. The resulting helio- 
seismic waves are observed as expanding circular ripples 
on the solar surface, which can be detected in Doppler- 
grams and as a characteristic ridge i n time-distance dia¬ 
grams dKosovichev fe Zharkovalll998l : IKosovichevI 2006l'l. or 
by ca l culating the integrated a coustic emission ( Donea et aid 
1 19991 : iDonea fe Lindsevl l2005l l . These flare-excited oscilla¬ 
tions are mostly local seismic waves. While the theory pre¬ 
dicts that global acoustic waves should also be excited, their 
amplitudes are thought to be significantly low er than the 
am plitudes of stochastically excited oscillations dKosovichevI 
l2009lb Most of the emitted acoustic energy of sunquakes is 
abo ve the acoustic cut- o ff freq uency of the Sun. 

iKaroff fe Kieldsenl d2008h concluded that global high- 
frequency solar acoustic waves have larger a mplitudes af¬ 
ter so me solar fl ares, a finding confirmed by iKumar et aid 
d20ldh However, I Richardson et ahl d2012[ ) found that a de¬ 
crease in acoustic power after a solar flare is just as likely 
as an increase. This is perhaps not surprising because the 
effect of the impulse on a global seismic mode depends on 
the location and time of the impulse. The effect of the im¬ 
pulse is therefore expected to increase the amplitude of some 
modes and decrease the amplitude in other modes. The re¬ 
sult should be a higher variation in amplitude distribution 
after a flare. 

In the case of the vastly more powerful stellar super¬ 
flares, the impact on the star may be much greater and 
it is conceivable that these forced oscill ations may b e ob- 
served. Thi s idea has been discussed bv iKaro 2 & and 
IKosovichevI (l2014l l in Kepler photometry of solar-like stars. 
Karoff found no significant enhancement of the energy in 
the post-flare acoustic spectra relative to the pre-flare en¬ 
ergy. However, a larger variability between the energy in the 
high-frequency part of the post- and pre-flare acoustic spec¬ 
tra was found compared to spectra taken at random times. 
This may be a result of the increased dispersion in acoustic 
energy discussed above. 

In spite of the very large differences in energy and op¬ 
tical emission between solar and stellar flares, the possibil¬ 
ity exists that the physical mechanisms for QPP could be 
very similar. In this paper we investigate QPP in Kepler su¬ 
perflares because this could provide additional clues to the 
nature of stellar flares. We also investigate the possible oc¬ 
currence of starquakes as a result of a superflare and search 
for excitation of global acoustic oscillation modes in the pe- 
riodogram. These investigations are only possible due to the 
superb precision of the Kepler data, a unique resource which 
is unlikely to be equaled for many years to come. 


2 DATA AND ANALYSIS TECHNIQUE 

Our analysis is based on th e 3140 flares in 290 stars observed 
in SC mode as described in IBalonal (120151 ). Out of these 3140 
flares, we selected 257 flares in 75 stars with sufficiently high 
signal-to-noise (S/N) to be suitable for detecting possible 
flare oscillations. 

In order to enhance the visibility of possible oscilla¬ 
tions it is necessary to remove the underlying flare decay. 
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Figure 1. Examples of four flares showing bumps on the decay branch. The top panel shows the flare light curve and fitted polynomial. 
The middle panel shows the decay branch with polynomial removed, with time measured from the time of maximum flare intensity. The 
bottom panel shows the wavelet spectrum as a function of time relative to the time of maxim um flare inten sity. The intensity is measured 
in parts per thousand (ppt). The KIC number and the flare identifier from the catalogue of lBalonal <|2015h is shown. 


We found that a polynomial of the form 

log (y) = do + ait + a 2 t 2 + 03 1 3 ... 

generally provides a good fit to the flare decay intensity. Here 
t is the time measured from flare maximum and y is the flare 
intensity. In nearly every case, a cubic polynomial was used. 
The polynomial was removed and the residuals plotted as a 
function of time. In about 18 percent of the Kepler flares, one 
or more bumps are present in the decay branch. A low-degree 
polynomial does not adequately remove these features. In 
these cases, the large-scale structures can be removed by 
applying a suitable filter to the data, as described below. 
Flare loop oscillations (QPP) are usually seen as 


damped sinusoidal oscillations (there are however decay¬ 
less, multi-modal and wave-train regimes observed). Wavelet 
analysis is an appropriate tool for detecting such a 
signal. Wavelet analysis decomposes a time series into 
time/frequency space simultaneously. One gets information 
on both the amplitude of any periodic signals within the se¬ 
ries, and how this amplitude varies with time. This is accom¬ 
plished by selecting an appropriate wavelet function that can 
be scaled and translated. The resulting wavelet transform is 
a representation of the signal at different time scales. We use 
the Morlet wavelet (a sine wave multiplied by a Gaussian en¬ 
velope) as the basis function in our analysis. The results are 
shown in a plot of instantaneous wave period and amplitude 
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Table 1. Physical parameters for stars in which one or more bumps are visible in the flar e light curve. The second column, Nf, is an 
identifier for the flare, which is the IVy-th flare for the particular star in the catalogue of iBalonal ("2015' ). The Kepler magnitude, Kp, 
effective temperature, T e g (K) and surface gravity, log g, is taken from the Kepler Input Catalogue (KIC). The luminosity, L/Lq is 
derived from the effective temperature and radius listed in the KIC. The average nu mber of flares per day, N/d, and the total number of 
observed flares, N, are shown. The rotation period, P TO t (d), are from lBalonal ll 201 ol i. P max is the expected period of solar-like oscillations 
(min). The last column is the derived period from wavelet analysis. 


KIC 

N f 

Kp 

(K) 

T e ff 

(dex) 

log g 

log L/Lq 

N/d 

N 

(d) 
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-0.9049 
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as a function of time. This is conveniently accomplished by applied to a data window with a length of a few days and 

a greyscale representation of the amplitude in a period-time starting a few hours after the flare. 

diagram. 

Forced global acoustic oscillations caused by a flare im¬ 
pulse are expected to be seen several hours after the flare 
and with a decay time of perhaps a few days. For detect¬ 
ing such oscillations we use standard periodogram analysis 
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Figure 2. The top panel shows the flare light curve and fitted 
spline curve. The second panel shows the decay branch with fitted 
curve removed. The third panel shows the wavelet spectrum. 
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Figure 3. The same as Fig.[2] 

3 FLARES WITH BUMPS IN THE LIGHT 
CURVE 

Of 257 stellar flares selected because of their very high 
S/N, 47 flares (i.e. 18 percent) show obvious structures (i.e. 
bumps) in the flare decay light curve. A fur ther 14 percen t 
show a distinct change in the rate of decay (lBalonall2015l) . 
Examples of flare light curves showing bumps are shown in 
Fig. [I] One possible reason for such bumps is that two flares 
occur by chance within the same small interval. One can 
estimate the chances of such an occurrence by measuring 
the mean flare rate during a suitable interval. The typical 
duration of a flare is about one hour. The star with the max¬ 
imum flare rate produces, on average, one flare every 10 hr 
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Figure 4. The same as Fig. [2] 


or x = 0.1 flare hr 1 . The probability, P, that two flares will 

occur in the interval of one hour follows a Poisson distri- 
2 _ 

bution and is P = ^-e x = 0.005. This is an upper limit 
considering the fact that we have chosen the maximum flare 
rate. In some stars (e.g. KIC 11551430) there are many flares 
with bumps. The probability of such an occurrence, assum¬ 
ing that what we are seeing is just a superposition of two 
or more practically simultaneous flares, is negligible. Hence 
there must be a physical process involved in producing the 
bumps. One such process could be sympa thetic flaring wher e 
one solar flare may trigger another flare (|Moon et al.ll2002l ). 

Both bumps and changes in rate of decay can be mod¬ 
eled as QPP with a rapid decay rate. In Fig. [T] the polyno¬ 
mial fit to the decaying branch, and the wavelet spectra are 
shown. Wavelet analysis of 47 “bump” flares in 30 stars in¬ 
dicates that significant power is present at a certain period. 
If we are to make progress in understanding this relatively 
common phenomenon, the “period” derived in this may offer 
some clue as to the nature of the bumps. 

If we assume that the bump is due to a highly-damped 
global acoustic mode, for example, one may expect the pe¬ 
riod to be correlated with some stellar parameter or combi¬ 
nation of parameters. Solar-like modes in stars are excited 
because their periods are similar to the typical turn-over 
period of a convective cell and below the acoustic cut-off 
frequency. As a result, the frequency of maximum ampli- 
tude, r ma x is related to the stellar paramet ers as follows 
(iBrown et al . 199ll : iKieldsen fe Beddinelll995l ) : 


^max ~ ^max© " 


M/Mq 


(R/R @ )WT e s/T em 
where the solar value for the frequency of maximum am¬ 
plitude is i/maxo = 3120/rHz, while M/Mq, R/Rq and 
Teff/T’efjQ is the stellar mass, radius and effective temper¬ 
ature relative to the Sun. 

If, by analogy with sunquakes, the stellar flares excite 
mostly acoustic modes with frequencies close to the acoustic 
cut-off frequency, one might expect a correlation between the 
period and the acoustic cut-off frequency. Since the acoustic 
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Figure 5. The same as Fig. [2] 




Figure 6. The same as Fig.[2] 


cut-off frequency is proportional to u ms . x , one may expect a 
correlation between the bump period and P max = l/i'max- 
The value of P max , as calculated from log g and T e g, is shown 
in Table[T] 

As shown in Table[2] the flare bump period (Pb um p) can 
vary widely even in the same star. We can find no correla¬ 
tion between this period and P max . In fact, no correlation 
between these timescales and any stellar parameter can be 
found. It seems that the structures on the decaying branch 
of the flare, whatever their cause, cannot be attributed to 
highly-damped impulsively excited global acoustic oscilla¬ 
tions. It is still possible that these may be highly-damped 
flare loop oscillations, but there is no evidence that can be 
used in support of this notion either. However, the similarity 
of the observed decaying QPP in stellar superflares, and of 


the standi ng oscillations observed in ho t coronal solar flare 
loops (e.g. IWangll201l| ; iKim et al .1120121 ) indicates the pos¬ 
sible similarity of the physical processes involved, despite 
differences in emission spectrum. 


4 RAPID OSCILLATIONS 

A polynomial fit to the flare decay branch is no longer ad¬ 
equate to enhance the visibility of possible oscillations with 
periods of only a few minutes. We found a segmented spline 
fit to be adequate for this purpose. In this method, a number 
of evenly spaced points are selected in the decay branch of 
the flare light curve. The mean intensity in the neighbour¬ 
hood of each point is found. A spline fit is calculated using 
the time and mean intensities at these points and the fit re¬ 
moved from the data. It is important to choose the interval 
between the segments carefully. If the interval is too small, 
possible oscillations with periods longer than this interval 
will not be detected because less than one period will be 
sampled. On the other hand, if the interval is too long, then 
the spline interpolation is no longer an adequate fit and the 
residuals may be contaminated by spurious long-period sig¬ 
nals. We found that a choice of 10-15 min was appropriate 
for the interpolation sampling interval in most cases. We 
have also tried other methods for trend removal, such as 
temporal smoothing, but this did not significantly influence 
the detected periods. 

We carefully examined all 257 flares without finding any 
obvious QPP in the majority of flares. We did not find any 
cases where QPP begins before the time of flare maximum. 
This may be due to the long sample time of 1 min and the 
steep rising branch which which would make such a detec¬ 
tion difficult. Evidence for damped oscillations after flare 
maximum is, however, present in the few flares shown in 
Figs-UHSl The stellar parameters for these stars are shown 
in Tablc[2] Most of these flares also appear in Tablc[l] other 
star s in this table do not show more than a single bump. 

iGruber et alj (1201 lh tested the apparent QPP of four 
bright solar flares observed in gamma rays using classical 
periodogram analysis, but found th at these os cillations were 
not intrinsic to the flares. Similarly. IVaughanl (120101 ) applied 
Bayesian statistics to apparent QPP in some Seyfert galax¬ 
ies and also found that these were not significant. The dif¬ 
ficulty is that the underlying noise in the periodogram is 
not ‘white’ (i.e., independent of frequency) but ‘red’ (i.e. in¬ 
creases towards low frequencies). This variation of noise with 
frequency needs to be taken into account in any analysis of 
significance. Because white noise is usually assumed the sig¬ 
nificance of QPP has generally bee n overestima t ed. T his is 
further illustrated by recent work bv llnglis et al.l (120151 ) who 
investigated supposed QPP in a selection of solar flares from 
a variety of sources as well as QPP in some optical stellar 
flares. They found that for all except one event tested, an 
explicit oscillation is not required in order to explain the 
observations. Instead, the flare signals are adequately de¬ 
scribed as a manifestation of a power law in the Fourier 
power spectrum, rather than a direct signature of oscillat¬ 
ing components or structures. 

That does not mean that all QPP signals are spurious, 
but that great caution needs to be exercised in determining 
the significance of QPP in flares. For this reason, we do 


© 0000 RAS, MNRAS 000, 000-000 

































Oscillations in stellar superflares 7 



Figure 7. The same as Fig.[2] 
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Figure 8. The same as Fig.[2] 


Table 2. Physical parameters for stars in which a damped oscilla¬ 
tion (QPP) is visible in the flare light curve. The second column, 
Nf, is an identifier for the flare. The Kepler magnitude, Kp, effec¬ 
tive temperature, T e fj (K) and surface gravity, log g, is taken from 
the Kepler Input Catalogue (KIC). The average number of flares 
per day, N/d, and the total number of observed flares, N, are 
shown. The variab ility type and the rotation period, P ro t (d), are 
from lBalonal < 12015 '). The last column is the period (min) derived 
from wavelet analysis. 
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solar flux tubes and solar flare loops. For a standing oscil¬ 
lation in a loop, the loop length L is given by L = jcP/2, 
where P is the oscillation period, j the parallel mode num¬ 
ber and c the appropriate wave speed. The waves with the 
longest periods, which are those of interest to us due to the 
1-min cadence of the Kepler data, are the slow modes. For 
slow modes, c is the tube speed, Ct with 



where c s is the sound speed and ca is the Alfven speed. 
Since the Alfven speed in the coronal loops is considerably 
larger than the sound speed, we assume c « c s . 

The unknown values of j, c and L render the extrac¬ 
tion of any meaningful physics impossible at this stage. We 
also note that the above description applies to a low-density 
plasma environment as occurs in a typical solar coronal flare 
loop. Stellar flares emit in a continuum and the physical 
process is likely to be different from that of solar flares even 
if the energy source, magnetic reconnection, is the same. 
However, these parameters are contained in the characteris¬ 
tic parameters of the decay phase of the flare, the damping 
time of the oscillations, intensity of the flare, and other ob¬ 
servables and, in principle, can be extracted from the data 
when a sufficiently detailed model becomes available. 


not claim that the oscillating signals seen in Figs.[2][8] are 
necessarily real. We merely wish to illustrate that QPP in 
the Kepler flares, if it exists, is not common. 

While QPP is a useful diagnostic tool for flares on the 
Sun, it cannot be used for stellar flares without some un¬ 
derlying assumptions. In the Sun one can image the flare, 
so that the loop length is known. Very often, other parame¬ 
ters, such as the plasma temperature, can also be quite well 
estimated. For stellar flares, however, we only have the pe¬ 
riod and decay time of the QPP, so the information that can 
be e xtracted is sev erely limited. IZaitsev fc Stepanov! (1 19821 ) 
and iRoberts et al.l dl984l ) showed that for a simple cylin¬ 
drical magnetic flux tube, several types of magneto-acoustic 
wave modes are possible: the slow (acoustic) mode, the fast 
kink and the fast sausage modes. These are all observed in 


5 FLARES IN STARS WITH SOLAR-LIKE 
OSCILLATIONS 

We know t hat some solar flares may exci t e local seismic dis - 
turbances dKosovichev fe Zharkova|[i998l : lKosovichevlf2006l ). 
with maximum power in waves of high spherical harmonic 
degree, l. These waves would not be visible in whole-disk 
photometry because the net brightness change is canceled 
owing to the small spatial wavelengths. The po wer in waves 
wit h low l which might be visible is very low dKosovichev! 
120051) . The more powerful impulse provided by a stellar 
superflare may, however, lead to observable results (star- 
quakes). As described in the Introduction, the expected re¬ 
sult is that the amplitudes of some global modes will grow 
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Figure 9. Periodograms of flare stars with solar-like oscillations. 


Table 3. Parameters for flare stars with solar-like oscillations. 
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while others will diminish after a flare. It is therefore to be 
expected that the net result will be a larger variability of the 
energy in the global modes after the flare which can manifest 
itself as a larger variability in the periodogram. 

The Kepler SC data is an ideal data set for studying 
possible starquakes. Among the 209 flare stars observed in 
SC mode, we can identify eight stars (Tablc[3J which clearly 
show the characteristic Gaussian envelope of solar-like os¬ 
cillations in the periodogram. The relevant portion of the 
peri odograms are shown in Fig. [9] 

iKaroffl (|2014|) calculated the pre- and post-flare acoustic 
spectra from substrings of different lengths before and after 
the flares. The photometric variability associated with the 
flares was evaluated by measuring the total energy in the 
high-frequency part of both the pre- and post-flare acous¬ 
tic spectra. Measuring the total power in a given frequency 
range may not be the most efficient method of detecting 
changes in the solar-like oscillations because the spectra are 
noise limited. In other words, a significant part of the total 
power comes from the noise and not the actual oscillations. 
Instead, it might be best to compare the amplitudes of in¬ 
dividual peaks in the periodogram of the global oscillations 
before and after the flare. In this way one may hope to de¬ 
tect possible systematic increases or decreases of amplitude 
in individual modes as might be expected. 

In Fig. [10] we show the light curve and periodogram of 
KIC 4554830 around the time of a flare. The periodogram is 
calculated using two or three 5-d data windows before and 
after the flare. Because of the stochastic nature of solar¬ 
like oscillations, changes in amplitude of individual peaks 
are to be expected. One therefore needs to examine several 
different frequency peaks and to determine if the changes in 
amplitude (increase or decrease) is significantly larger than 
would normally occur. In KIC 4554830 the amplitude of the 
mode at 1790.99 /rHz appears to increase after the flare. On 
the other hand, the amplitude of the mode at 1839.41 /j ,Hz 
seems to decrease. The data windows are independent (no 
overlap) so one can judge the significance of the amplitude 
changes. 

Similar diagrams may be created for all flare stars with 
solar-like oscillations, but this may not be the best method 
of detecting starquakes. Since we are looking for changes in 
amplitude of individual modes, we selected several peaks of 
large amplitude and fitted a truncated Fourier series to the 
data window using these frequencies. This gives us the am- 
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Figure 10. Top panel: light curve of KIC 4554830 showing flare 
and regions where periodograms were calculated. Bottom panel: 
periodograms of the regions shown in the top panel showing the 
solar-like oscillations. The dotted line shows the time of the flare. 


plitude and its standard deviation for each frequency peak. 
We performed this calculation in two data windows before 
and after the flare. 

Fig-ED shows the resulting amplitudes in each 5-d win¬ 
dow as a function of time for all solar-like modes of sufficient 
amplitude. For clarity, all amplitudes are normalized to their 
values in the window just before the flare. The same mode is 
connected with a solid line. Judging from the error bars, it 
is evident that in no case is there a clear amplitude change 
after the flare relative to the pre-flare amplitude. We con¬ 
clude that a superflare has no systematic influence on the 
amplitudes of the solar-like oscillations detectable by this 
technique. 

This does not mean that starquakes do not occur, of 
course. The impulse generated by a superflare is bound to 
cause acoustic disturbances which could affect the ampli¬ 
tudes of global modes. However, the small spatial scale (high 
azimuthal degree) of such oscillations lead to cancellation 
effects, resulting in whole-disk light variations that are too 
small to detect. 
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6 FLARES IN <5 SCUTI AND 7 DORADUS 
STARS 

It is very difficult to detect flares in classical pulsating 
stars because the light variations due to the pulsation tend 
to mask short-lived rapid excursions such as a flare. This 
is especially true for short-period variables. Nevertheless, 
among the SC observations, flares are seen in the d Set stars 


Figure 11. Pulsation amplitudes for solar-like oscillations as a 
function of time. The amplitudes have been normalized to the 
values just before the flare and one-er error bars are shown. Solid 
lines join individual mode frequencies. The time of the flare is 
indicated by the vertical line. 
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Figure 12. Periodograms of a 7 Dor star (KIC 5113797, top 
panel) and two 5 Set stars that show flares. 


KIC 1294756 and KIC 2301163. Both these stars have ex¬ 
tremely low pulsation amplitudes, rendering the flares more 
easily visible. The flare star KIC 5113797 seems to be a 7 Dor 
variable of low amplitude. The periodograms of these stars 
are shown in Fig. 1121 

The oscillations in S Set and 7 Dor stars are self-excited 
with reasonably stable amplitudes and phases, unlike the 
stochastic solar-like oscillations. Hence changes in pulsation 
amplitude after a flare might be easier to detect. Unfortu¬ 
nately, the S/N ratio in the 6 Set pulsations of KIC 2301163 
is too low for the oscillations to be detected in a 5-d data 
window. Since the effect of a flare impulse on the oscillations 
is expected to dissipate rather quickly, one needs to use a 
short data window to optimize its detection. For this reason 
KIC 2301163 was excluded from the analysis. The low oscil¬ 
lation frequencies in the 7 Dor variable KIC 5113797 means 
that very few pulsation cycles can be obtained during the 
5-d window. As a result, the amplitudes and phases have 
large errors. This star, too, was omitted. 

In Fig.d part of the light curve of the $ Set star 
KIC 1294756 is shown centered on one of the flares. The 
boxes show the data windows used to construct the peri¬ 
odograms. Some amplitude changes seem to occur, particu¬ 
larly for the mode at 34.9317 d -1 which appears to decrease 
after the flare. 

In Fig.[l4]we show how the pulsation amplitudes vary 
with time before and after the flare. In this figure the data 
for two flares are shown. As before, the amplitudes are rela¬ 
tive to the amplitudes just before the flare. Although there 
are variations in amplitude for many modes, these are within 
the expected errors and are therefore not significant. We con¬ 
clude that even when the pulsation amplitudes are stable, 
there is no evidence that a flare affects the mode amplitudes 
in a manner that is detectable. 
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Figure 13. Top panel: light curve of the <5 Set star KIC 1294756 
showing flare and regions where periodograms were calculated. 
Bottom panel: periodograms of the regions shown in the top panel 
showing the 5 Set pulsations. 
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7 CONCLUSION 

In this paper we investigate possible periodic light variations 
arising from stellar flares. By analogy to the Sun, quasi- 
periodic pulsations (QPP) may be a result of MHD forces 
and other processes operating in the flare loop. These oscilla- 


Figure 14. Relative pulsation amplitudes for the <5 Set star 
KIC 1294756 as a function of time. The amplitudes have been 
normalized to the values just before the flare and one-cr error 
bars are shown. The time of the flare is indicated by the vertical 
line. 
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tions, which are often observed in solar flares, have also been 
detected in other stars and have periods in the range of sec¬ 
onds to tens of minutes. The period of the oscillation offers 
the potential to probe the magnetic field strength and tem¬ 
perature in the flare. However, nearly all the light in solar 
QPP is emitted in lines of highly-ionized elements, whereas 
observations of stellar QPP are essentially white light. One 
cannot therefore directly compare solar and stellar QPP, 
though the underlying mechanism might still be the same. 

QPP in some stellar flares obser v ed from the ground 
(|Rc^onj_[l974|^_|ZM alJ_[2000||__ [Mathio^ et al.l 


l2006l : Contadakis et al. 201dl : Qian et al.l 2012l l have been 
interpreted as analogous to flare QPP in the Sun, in spite 
of the very different light emission properties described 
above. These QPP generally have very short periods and 
last for several cycles. We searched for short-period QPP 
using wavelet analysis in all flares with high S/N and found 
seven flares in five stars in which such an oscillation seems 
to be present (Figs.[2][8]). The periods do not correlate with 
v max or any stellar parameter. These are probably examples 
of the same phenomenon seen in the ground-based observa¬ 
tions discussed above. They may not be directly comparable 
with QPP in solar flares, but perhaps a similar process may 
be active. 

Oscillations in the light curve could also arise as a result 
of the impulse generated by a flare (a starquake). In the Sun, 
these are seen as surface waves of short spatial wavelength 
radiating from the location of the flare. There is evidence 
that the impact also affects the amplitudes of the stochastic 
global oscillations. For some modes a starquake will lead to 
an increase in amplitude, while for other modes a decrease 
in amplitude may be expected. 

A significant fraction of stellar flares show clear struc¬ 
tures (bumps) on the decaying branch of the light curve. 
We argue on probability grounds that these cannot be a re¬ 
sult of almost simultaneous multiple flares. The bump may 
be modeled as rapidly decaying QPP. We tested the possi¬ 
bility that the bumps may be highly damped forced global 
oscillations by measuring the period using wavelet analysis. 
If these are forced global oscillations, one might expect the 
period to correlate with some stellar parameter. We could 
find no correlation with the estimated frequency of maxi¬ 
mum amplitude, i^ max , resulting from solar-like oscillations 
or any other stellar parameter. We conclude that the bump 
cannot be understood as a forced global oscillation. 

Finally, we attempted to detect forced global oscilla¬ 
tions resulting from a starquake. These are expected to be 
visible shortly after a flare and are expected to modify the 
amplitudes of individual modes in stars where solar-like os¬ 
cillations are detected. The best way to look for this ef¬ 
fect is to measure the amplitudes of individual modes be¬ 
fore and after a flare. There are eight stars in the Kepler 
short-cadence observations which show flares and solar-like 
oscillations. Examination of the amplitudes of modes before 
and after a flare showed no obvious indication of significant 
amplitude changes. We conclude that the effect is too small 
to be detected in the Kepler data. 

By their nature, random changes in amplitude are a 
characteristic of solar-like oscillations and this may mask an 
amplitude change resulting from a flare. The self-excited os¬ 
cillations in 5 Scuti and 7 Doradus stars lead to generally 
stable amplitudes. Therefore these stars may provide better 


detection of starquakes. We analyzed the pulsation ampli¬ 
tudes of two 5 Set and one 7 Dor star before and after a 
flare. Again, we were not able to find evidence of amplitude 
variations due to a flare. 

A major problem in our understanding of stellar flares 
is that there are no corresponding observations in the Sun. 
Apart from the huge disparity in energy, solar flares emit 
almost entirely in emission lines of highly-ionized elements 
whereas stellar flares are essentially white light flares. Fur¬ 
thermore, there are no reported observations of QPP in 
white-light solar flares. Such Q PP may be detected i n 
stellar flares observed in X-rays llMitr a-Kr aev et al.,1120051 ). 
The standard flare model suggests that the white-light 
emission in solar and stellar flares is triggered by non- 
thermal electrons which originate in the corona. There is 
a strong corresponde nce between white-l ight emission and 
hard X-ray emission (iHudson et al.ll2006h . Moreover, solar 
hard X-ray bursts often show a high degree of periodicity 
(lAschwanden et ai1ll994 ). It is possible that the QPP seen 
in the few Kepler flares may be a result of these beams and 
their effect in the lower atmosphere. However, it may be dif¬ 
ficult to understand white light flares in A stars where it is 
generally assumed that a corona is not present. 

Our conclusion is that white-light QPP may possibly 
be seen in some Kepler flare stars, but their nature may dif¬ 
fer from QPP in solar flares, although the processes involved 
could be similar. It would be important to synthesize whole- 
disk white light observations of solar flares. In this way we 
may hope to extend what we know of solar flares to stel¬ 
lar flares and thereby create a fuller understanding of the 
mechanisms involved in solar and stellar flares. 
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